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We present the results of an ab initio study of magnetic properties of Fe, Co and Ni surfaces. 

In particular, we discuss their electronic structure and magnetic exchange interactions (Jij), as 
obtained by means of a combination of density functional theory and dynamical mean-field theory. 

All studied systems have a pronounced tendency to ferromagnetism both for bulk and surface atoms. 

The presence of narrow-band surface states is shown to enhance the magnetic moment as well as 
the exchange couplings. The most interesting results were obtained for the Fe surface where the 
atoms have a tendency to couple antiferromagnetically with each other. This interaction is relatively 
small, when compared to interlayer ferromagnetic interaction, and strongly depends on the lattice 
parameter. Local correlation effects are shown to lead to strong changes of the overall shape of 
the spectral functions. However, they seem to not play a decisive role on the overall picture of the 
magnetic couplings studied here. We have also investigated the influence of correlations on the spin 
and orbital moments of the bulk-like and surface atoms. We found that dynamical correlations in 
general lead to enhanced values of the orbital moment. 

PACS numbers: 75.30.Et, 31.15.E-, 71.27.+a, 73.20.At 


I. INTRODUCTION 

Bulk Fe, Co and Ni are all classical examples of ferro- 
magnets. However, when confined to two dimensions, 
these transition metals (TM) show a large panorama 
of fascinating magnetic properties and phenomena [1 
3]. For instance, thin layers of these atoms may show 
antiferromagnetic (AFM) behavior or even non-collinear 
spin structures, depending on the film thickness and/or 
the substrate [4-6]. The latter is known to play an im¬ 
portant role, producing strains due to lattice mismatch 
and hybridizing with the TM states [7-10]. All above- 
mentioned effects contribute to the magnetic exchange in¬ 
teractions (Jij), that are the relevant parameters of an ef¬ 
fective spin-Hamiltonian which determine the Curie tem¬ 
perature and magnon dispersion of the material. The lat¬ 
ter two quantities are of particular importance for tech¬ 
nological applications in, e. g., spintronic memory and 
logic devices. Hence, a fundamental understanding of the 
magnetic properties of these systems is needed. 

The magnetism of surfaces has been of interest for quite 
some time. Initial studies were mainly focused on differ¬ 
ences between surfaces and bulk properties. For instance, 
experimentally it was for some time discussed that fee Ni 
layers on top of a Cu substrate produced magnetically 
’dead’ layers, with an absence of magnetic moments [11]. 
However subsequent experiments [12] and theory [13-15] 
suggested that the spin-moments at surfaces in general 
are enhanced, since the bands are narrower. Later on rel¬ 
ativistic electronic structure theory could analyze also the 


orbital moments of surfaces, and here the enhancement of 
the surface magnetism was found to be even larger than 
the spin-contribution for bcc Fe, hep Co and fee Ni [16- 
18]. These theoretical predictions were confirmed by ex¬ 
periments using x-ray magnetic circular dichroism [19]. 

Computational modeling is very important to inves¬ 
tigate magnetic properties, as it gives a material spe¬ 
cific description and makes it possible to disentangle all 
the relevant contributions. Density functional theory 
(DFT) and its formal extensions [20-22] give an excel¬ 
lent parameter-free description of ground-state proper¬ 
ties of magnetic metals, including bulk structures as well 
as systems without three-dimensional periodicity such 
as surfaces, interfaces, thin films, disordered alloys and 
nanoparticles. However, several studies have emphasized 
the importance of including strong correlation effects in 
the electronic structure of bulk Fe, Co, and Ni. For in¬ 
stance, noncoherent features such as Hubbard bands and 
satellites which appear in the excitation spectra of the 
photoemission experiments [23, 24] can not be described 
by LDA and/or GGA. In addition, LDA calculations pre¬ 
dicts too wide majority spin 3 d band and overestimate 
the spin splitting for these materials [25-28]. 

Correlation effects in transition metals are expected 
to be even more pronounced for the surface atoms, due 
to narrower bands and reduced coordination numbers. 
In this article we report on a computational study of 
surface magnetism of TM slabs. The main focus is on 
the calculations of magnetic moments and inter-atomic 
exchange interactions (Jij)- The simulations are based 
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on a combination of DFT and dynamical mean-field the¬ 
ory (DMFT). This technique, which is usually addressed 
as LDA+DMFT, is at the moment the state-of-the-art 
method to study strong correlations in materials at finite 
temperature [29-31]. 

To the best of our knowledge, LDA+DMFT has not 
been previously applied to the exchange interactions of 
the transition metals surfaces. Even for standard DFT 
such simulations are rare, as most attention was focused 
on thin films on various substrates [10, 32-35]. One of 
the reasons for this is that many softwares are still based 
on the atomic-sphere approximation (ASA), which limits 
their use for studying surfaces, or low dimensional sys¬ 
tems in general. The methods used in the present work 
do not suffer of this limitation. 

This paper is organized as follows. In section II we 
briefly explain the computational scheme used in this 
work, as well as the implementation of the formalism 
described in Ref [36] for evaluating the exchange param¬ 
eters. The result of our electronic structure calculations 
and exchange interactions for Fe, Co and Ni surfaces are 
presented in section III. The following section reports our 
investigation of the orbital polarizations for each slab. Fi¬ 
nally, we draw our conclusions which will be followed by 
three Appendices. Appendix A concerns the influence 
of the Hubbard U value, Appendix B is about the ef¬ 
fect of full self-consistency over the charge density on the 
exchange parameters and Appendix C describes renor¬ 
malization factors due to many-body effects. 


II. THEORY 

The electronic structure as well as the magnetic prop¬ 
erties of the TM slabs were investigated in the frame¬ 
work of scalar-relativistic full-potential linear muffin-tin 
orbital (FP-LMTO) code RSPt [37]. Due to the full 
potential character the code does not have limitations 
dictated by the geometry of the problem under consid¬ 
eration. Moreover, due to the small number of basis 
functions, RSPt is particularly suitable for LDA+DMFT 
simulations with full self-consistence over self-energy and 
electron density. Details of this implementation were pre¬ 
sented elsewhere [27, 37 39] and will not be repeated 
here. We redirect the reader to those references for a 
detailed overview of our formalism. 

Once the electronic structure was converged, the 
magnetic excitations were mapped onto the Heisenberg 
Hamiltonian: 


H— 'y [ Jij e t Cj , (1) 

i+J 

where J- t j is an exchange interaction between the spins, 
located at sites i and j , and e) is a unity vector along the 
magnetization direction at the corresponding site. We 
extracted the pair-wise exchange interactions by employ¬ 
ing the method of infinitesimal rotation of the spins. The 


exchange parameters were computed using the local mag¬ 
netic force approach [36, 40], which reads: 


J — -Tr 

'Jij — ^ I u. 


X!(itUn)GUiu n )^(lCU n )GUiLOn) , ( 2 ) 


where Gfi is the inter-site Green’s function, a denotes 
spin projection (<r = {f,i}), and the trace is taken over 
the fermionic Matsubara frequencies iu> n and the states 
characterised by an angular momentum quantum number 
L. The crucial quantity in Eq. 2 is the dynamical on-site 
exchange potential: 

±t(ico n ) = ( H J - Hi) + (eJ (iu n ) - Et(iw n )) , (3) 


where H? is the local Hamiltonian matrix, obtained by 
solving the DFT equations and Ef is the self-energy, de¬ 
scribing the electronic correlations. The self-energy ap¬ 
pears only for LDA+DMFT calculations, and also enters 
the expression of the Green’s function as: 


Gij{iu) n ) = ( i 


1 


iu n - H - E (iui n ) 


(4) 


More details about the evaluation of the exchange inter¬ 
actions, in particular in relation to the basis set used for 
the local orbitals, can be found in Ref. 41. 

Relativistic effects will not be considered in our work, 
for sake of simplicity and unless explicitly stated. These 
effects give rise to other types of magnetic interactions, 
like anisotropic exchange couplings and magnetocrys¬ 
talline anisotropy. However, the bilinear term, consid¬ 
ered in the present work, is usually the leading one. 
For instance, it was recently shown for Fe/Rh(001) that 
by considering Heisenberg interactions only it is possi¬ 
ble to obtain a very detailed picture of magnetic excita¬ 
tions, giving an excellent agreement with experiment [34]. 
Nevertheless, we performed a few additional simulations 
with spin-orbit coupling included in order to analyze the 
enhancement of the orbital magnetism at the surfaces, 
which is a very important problem in materials science. 
These results will be presented at the end of the paper. 


Computational Details 

DFT simulations were performed by using LDA as 
exchange-correlation functional. After the convergence, 
we have applied a LDA+DMFT technique for a selected 
set of TM 3d orbitals. The k integration over the irre¬ 
ducible wedge of the Brillouin zone have been performed 
using 24x24x24 points for bulk and 24x24x1 points for 
the slabs. We have performed relaxation of the topmost 
layers of Fe slab, which are known to be quite small in 
these systems [42]. In LDA, for Fe, we obtained a 0.1% 
(1%) reduction of the surface (subsurface) magnetic mo¬ 
ments with respect to unrelaxed slabs (truncated bulk). 
For Co and Ni, we did not perform extensive tests as 
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we expect the changes induced by the relaxation to be 
even smaller. In fact these changes are proportional to 
the difference between bulk and surface spin moments, 
which is much larger for Fe than for Co and Ni. There¬ 
fore, to avoid presenting two sets of similar results and 
to facilitate comparison with similar studies, our analysis 
will be limited to unrelaxed slabs, where the interatomic 
distances depend solely on the bulk lattice parameter. 
The latter was chosen as obtained from experiments, i.e. 
2.86 A for bcc Fe, 3.52 A for fee Ni and 2.51 A for hep 
Co [43]. For the latter, the distance between the hexag¬ 
onal planes was chosen as 4.07 A [43]. The free-standing 
slabs of Fe, Co and Ni have been modeled using 15 layers 
of their bulk structure repeated in (001) direction for Fe 
and Ni, while the (0001) direction was used in the case 
of Co. Since three dimensional periodic boundary condi¬ 
tions are used, a 27-A-thick layer of vacuum was used to 
construct a supercell. 

LDA+DMFT simulations were performed for a tem¬ 
perature of 400 K. The effective impurity problem arising 
in DMFT was solved through the spin-polarized T-matrix 
fluctuation-exchange (SPTF) solver [44]. Since the lat¬ 
ter is a perturbative approach, it can only be applied to 
systems with moderate correlations and in the metallic 
regime of the Mott-Hubbard transition. SPTF is usually 
applied by using the static part of the self-energy as a 
double-counting correction term. This choice has been 
used for all DMFT simulations, throughout the paper. 

For 3d orbitals, where the electrons are supposed to 
show more atomic-like features, the Coulomb interaction 
can be parameterised via Slater integrals F rl [45]: 


U = F °, 


J = 


F 2 + F 4 
14 


(5) 


where U is Hubbard parameter and J is Hund’s exchange. 
The values of U and J can be either extracted from exper¬ 
iments or calculated from first principles. In this work, 
we have taken their values from the literature [25, 26, 44], 
LDA+DMFT calculations for Fe and Co were done uti¬ 
lizing U = 2.3 eV and J = 0.9 eV, while for Ni U = 3 eV 
and J = 0.9 eV were chosen. In order to see the effect 
of U on the spectra and exchange parameters, we have 
performed some test calculations using larger values of 
U for surface and subsurface atoms, while keeping the 
previous value unchanged for the inner layers. Although 
it is known that the choice of U affects the intensity and 
the position of satellites in the valence band spectra, e.g. 
as reported in Ref. [27], our results suggest that the ex¬ 
change parameters are marginally affected by varying the 
U value. This analysis is illustrated in Appendix A. 

Finally, we performed extensive tests to analyze the 
role of self-consistence over the electron density in the 
LDA+DMFT cycle. We focused on spectral functions 
and exchange interactions for 7-layer slab. These results 
are presented in Appendix B. Our general conclusion is 
that updating the electron density in the LDA+DMFT 
cycle introduces minor corrections for the transition met¬ 
als surfaces, at least concerning the magnetic properties. 



FIG. 1. Layer-resolved Fe-3d projected spin moment for dif¬ 
ferent slabs, whose thickness varies from 9 layers to 31 layers. 
The layers are numbered from 1 (surface) to A max (the inner¬ 
most layer). A max is equal to the total number of layers in 
the slab plus one, divided by two. 


In light of these minor changes as well as for compu¬ 
tational efficiency, we performed simulations of 15-layer 
slabs by keeping the electron density unchanged. 


III. RESULTS 


We consider free-standing slabs of 3 d transition metals 
in their most stable magnetic structure below the Curie 
temperature. The exchange interaction is calculated for 
different layers as a function of the interatomic distance 
between pair of atoms. These calculations are aimed at 
understanding in detail the differences between bulk and 
surface, and also at seeing how the local dynamical corre¬ 
lations affect the exchange interactions. In the next few 
sections we will elaborate on these issues separately for 
each element. 


A. Fe 

It is important to address first the convergence of the 
relevant magnetic properties with respect to the thick¬ 
ness of the slab. We present this analysis only for Fe for 
brevity. In Fig. 1 layer-resolved Fe-3d projected spin mo¬ 
ments are reported for slabs of different thickness. These 
calculations reveal long range damped oscillations of the 
local moments when going from the outermost layer to 
the innermost one. This behavior is due to the surface 
induced changes in the magnetism of the itinerant ferro- 
magnets (Friedel oscillations) [46, 47]. In principle, the 
formation of the quantum states in the finite size slabs 
is accompanied by the creation of a barrier on the sur¬ 
face which leads to a different electronic structure around 
the Fermi level (Ep). In the cases of Fe, Co and Ni, the 
magnetic moment of the innermost layer reached the bulk 
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FIG. 2. Layer-resolved projected density of states of 3 d or¬ 
bitals of the innermost layer (middle) of Fe slab and of the 
bulk for majority and minority spin components in LDA (left 
panel) and LDA+DMFT (right panel). 


value, with a difference smaller than 1%, for slabs of 15 
layers. For this size, the exchange interactions of the in¬ 
ner most layer of Co and Ni slabs, were equal to those of 
their bulk up to 0.5%. For Fe, however, a slightly larger 
difference was obtained, of about 1.2%, which was due 
to the difficulties in matching the same special points for 
the sampling of two-dimensional and three-dimensional 
Brillouin zones. Finally, an analogous convergence of the 
electronic structure can be observed in Fig. 2, where the 
projected density of states (PDOS) of the innermost layer 
is compared with the PDOS of the bulk for the 15-layer 
slab. The curves do not exhibit any visible difference on 
the scale of interest, for both LDA and LDA+DMFT. 

Next we have analyzed the differences in the PDOS of 
the surface atoms and that of the innermost layer. These 
results obtained in LDA and LDA+DMFT are shown 
in Fig. 3. The PDOS at the surface is very different 
from the bulk, due to a reduced coordination number, 
which results in narrower bands and more pronounced 
correlation effects. Our results are in a good agreement 
with prior studies reported by Grechnev et al. [27] and 
Chuang et al. [33]. Note that in the latter work the cal¬ 
culations were performed using generalized gradient ap¬ 
proximation, which accounts for some differences with 
respect to the results reported here. 

Fig. 3 shows that there are a large number of ma- 



FIG. 3. Layer-resolved projected density of states of 3 d or¬ 
bitals of Fe slab for the atoms at the surface, subsurface 
and the innermost layer for majority and minority spin com¬ 
ponents in LDA (left panel) and in DMFT approach (right 
panel). 


jority spin states in the vicinity of Ep for the inner¬ 
most layer of Fe slab, for both LDA and LDA+DMFT. 
These states arise mainly due to d yzi d xz , d xy orbitals (not 
shown here), and are shifted to lower energies for the 
atoms sitting at the surface. This results into an effec¬ 
tive suppression of spectral weight at Ep. For minority 
spin, instead, the bulk (innermost layer) PDOS shows 
just a few states of the t 2 g character around Ep , where a 
pseudo-gap forms. At the surface, the increase of the ex¬ 
change splitting causes these states to move just across 
Ep, which results in a drastic increase of the spectral 
weight. As a result, the spectral weight at the Fermi 
level arises mainly from one spin channel, which makes 
the surface behave as a strong ferromagnet. The inner 
layers have instead the characteristics of a weak ferro¬ 
magnet, similarly to the bulk. From Fig. 3, we also notice 
that the PDOS for the atoms sitting on the subsurface 
layer does not show substantial differences with respect 
to the PDOS of the bulk. 

We can now focus on the comparison between LDA 
and LDA+DMFT. Although the overall PDOS obtained 
by means of these two methods are quite different, they 
exhibit very similar behavior in the vicinity of Ep. Given 
that this region is of primary importance for the exchange 
interaction, we expect to obtain similar results within 
these two approaches, at least as concerns the asymptotic 
behavior. 

From this point on, we focus exclusively on the 15-layer 
slab. In Fig. 4, the layer resolved exchange parameters 
(Jij) are reported, for both LDA and LDA+DMFT. For 
clarity, we report only results for the most physically in¬ 
teresting layers, e.g. surface, subsurface and the inner¬ 
most layer. We will anyway make general considerations 
regarding all layers in the following discussion. The in¬ 
tralayer exchange interaction is referred to the case when 
the two atoms interacting with each other are located in 
the same layer. The interlayer interaction is referred to 
the case when the two atoms belong to different layers. 
The layers in the plots are denoted by 1 for the surface, 
2 for the subsurface and so on, analogously to Fig. 1. 

The first general consideration to draw from our calcu¬ 
lations is that interlayer exchange parameters for atoms 
in the inner layers are substantially smaller than those for 
atoms in layers closer to the surface. This trend is ob¬ 
served for both LDA and in LDA+DMFT. For instance, 
Fig. 4 shows that the exchange interaction between an 
atom at the surface and its first nearest neighbor (NN) 
sitting in the subsurface (blue lines in the right panel) 
is strongly ferromagnetic. The strength of this exchange 
interaction is about twice larger than that of an atom in 
the innermost layer and its first NN in an adjacent layer 
(pink lines in the left panel of Fig. 4). 

The quantitative explanation for the layer dependence 
of the exchange interactions can be provided by the two 
important factors. First, there is a direct influence of the 
on-site exchange field at the atoms i and j on the cor¬ 
responding parameter (see Eq. (2)). The second im¬ 
portant factor is the coordination numbers which would 
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FIG. 4. Layer-resolved exchange parameters (Jij) for 15-layer bcc Fe(001) slab for the case when atom i is located in the 
innermost layer (left panel), in the subsurface (middle panel) and at the surface (right panel). The solid lines indicate LDA 
results while the dashed lines represent the LDA+DMFT results. The layer numbering in the legend starts from surface denoted 
by 1, subsurface denoted by 2 and so on. The innermost layer is denoted by 8. The interaction between the surface and layer 
9 and further are small and are not shown in these plots. 


affect the Hamiltonian H , the self-energy and, therefore, 
the inter-site Green’s function (see Eq. (3)). 

An inspection of Fig. 3 has already revealed that ma¬ 
jority and minority spin states are more split for the sur¬ 
face atoms. Hence, the local exchange field (E s ) is overall 
larger than its bulk counterpart, which can explains the 
enhancement of the exchange integrals between surface 
and subsurface atoms with respect to couplings between 
more internal (adjacent) layers. Moreover, the surface 
PDOS is characterized by a larger spin polarization of 
the Fermi surface in comparison to inner layers. This 
provides a large number of available states with a certain 
spin projection right above Ep. Hence, similarly to the 
double exchange mechanism, an electron hopping will be 
facilitated if it does not have to flip its spin, i.e. if the 
neighboring moments are parallel to each other. This 
scenario thus also supports an enhancement of the ferro¬ 
magnetic interaction with the surface atoms, as obtained 
in our calculations. 

We then proceed to the analysis of the intralayer cou¬ 
pling. An interesting finding of the present investigation 
is that two Fe atoms at the surface possess an AFM cou¬ 
pling (green lines in the right panel of Fig 4). Note that 
this neighborhood corresponds to second NN atoms in 
bcc structure, and is therefore always smaller than the 
leading FM contribution between first NN. The presence 
of AFM exchange interactions at the surface is rather 
surprising, albeit we note that previous first principles 
theory also suggest such a coupling [32]. Fe is known to 
have the tendency to AFM coupling for hep or fee crystal 
structures [5, 7, 10, 48] as well as for thin mono-layers on 
some substrates, but seldom in bcc like environments [33- 
35]. For example, the exchange coupling between the two 
neighboring atoms at the surface of Fe clusters was re¬ 
ported to be FM [49]. 

Here we focus on the understanding the origin of this 


tendency to AFM coupling at the surface. For this pur¬ 
pose we analyze individual orbital contributions to the 
exchange parameter. The local Hamiltonian for each Fe 
atom is diagonal in the basis of cubic harmonics and so is 
E s from Eq. 2. Having E s in a diagonal form allows us to 
write each exchange coupling as J 12 = m2 J™ 1 ,rr ‘ 2 , 
where orbital m 1 is located at the site 1 and orbital m 2 is 
at the site 2. Exchange interaction between two closest 
surface spins in a form of a matrix in orbitals space is 
shown in Table I. The table hence shows the strength 
of the exchange interactions of symmetry resolved states 
of one atom with symmetry resolved states of a nearest 
neighbor surface atom. The total interaction between 
these two atoms is obtained by summing all components 
of Table I. The analysis reveals that there are basically 
two competing contributions to the between the NN’s 
surface moments. A first FM contribution originates 
from d xy — d xy bonds and d yz — d yz bonds, depending on 
the bond vector. A second AFM contribution, instead, 
arises from d yz —d x 2_ y 2 bonds. This contribution is much 
stronger in our case, and overcomes the FM part of the 
exchange. Similar competition was shown to take place 
for the next NN exchange couplings in bulk Fe, which cor¬ 
responds to the same coordination shell as we investigate 
here. [41] However, the balance between the two contri¬ 
butions can be easily changed by small changes of the 
NN distance. For example, Chuang et al. [33] reported a 
FM coupling between the adjacent surface atoms in free 
standing Fe films when using the Ir lattice constant. We 
have repeated their calculations and obtained the same 
results. We conclude that the tendency to AFM coupling 
is innate in the magnetic properties of the Fe surface. 
This suggests that the primary role of the substrate con¬ 
sists in modifying inter-atonric distances and not in af¬ 
fecting the electronic structure via direct hybridization, 
at least for the aforementioned cases. 
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FIG. 5. Panels (a) and (b): the calculated band structure 
of Fe slab for majority and minority-spin states, respectively, 
shown with the amount of 3d orbital character of the surface 
atoms (red color). Panel (c): The Fermi surface cross-section 
of the minority-spin bands, having the largest contributions 
from surface 3d states. Suggested nesting vector is indicated 
with the black arrow (the spectra are calculated in the com¬ 
plex energy E + i5 with 5=0.005 Ry). 

In order to have a more clearer explanation of the AFM 
interactions at the Fe surface, we have studied its under¬ 
lying electronic structure in more detail. In Fig. 5(a,b) we 
show the spin-polarized band structure, projected onto 
3d states of the surface atoms. The largest contribu¬ 
tion of these orbitals to each band is shown with the red 
colour. As one can see from Fig. 5, the surface states con¬ 
tribute only to the majority-spin bands well below Ep. 
In contrast, there is a significant weight of them at Ep 
for the minority-spin channel, in particular near T and A' 
points, which will contribute to the magnetic susceptibil¬ 
ity. To strengthen this point, in Fig. 5(c) we show parts 
of the Fermi surface cross-section, originating from the 
bands with a large surface component. One can see that 
the Fermi surface is nested between T and X points and 
the nesting vector is indicated with an arrow in Fig.5(c). 
The nesting vector (indicated with the black arrow) is di¬ 
rected along (100) direction and has a length very close to 
7r/a. This defines a preferable direction for a symmetry¬ 
breaking in the system and we suggest it is connected to 
the AFM coupling between the NN spins at the surface. 
We note as well that bulk bcc Fe has similar features of 
the Fermi surface (see e.g. Fig. 3 in Ref. 50). However, 
the nesting vector connecting the pockets located at T 
and H points has much a smaller length, which does not 
lead to the pronounced AFM interaction. 


TABLE I. Orbital-decomposed exchange interaction parame¬ 
ter between two closest neighbors at the surface for the bond 
vector (010). The reported values are in meV. 



d z 2 

d x 2 —y 2 

dyz 

dxz 

dxy 

d z 2 

0.81 

-1.02 

0.36 

0.00 

0.00 

d x 2_ y 2 

-1.02 

1.26 

-6.26 

0.00 

0.00 

dy Z 

0.36 

-6.26 

3.99 

0.00 

0.00 

dxz 

0.00 

0.00 

0.00 

-1.07 

0.05 

dxy 

0.00 

0.00 

0.00 

0.05 

6.32 


Other considerations can be drawn from the calculated 
exchange parameters from Fig. 4. AFM intralayer cou¬ 
pling between the closest neighboring atoms at the sur¬ 
face (light-green lines in the right panel) might in prin¬ 
ciple lead to non-collinear spin configurations. However, 
we estimate the resulting frustration to be weak, since 
their interlayer couplings with atoms in the substrate are 
much larger (dark-blue lines in the right panel). 

At the surface, the AFM exchange coupling between 
the closest neighbors as well as the lower coordination 
number result in smaller values for the total exchange 
interactions J\ = JT J\j . This consequence is, quali¬ 
tatively, in agreement with the LMTO-ASA results re¬ 
ported by Turek ef. al [32]. A quantitative comparison 
is, however, not possible for Fe, for which the long-range 
oscillatory (RKKY-like) behavior of exchange constants 
makes the value of J\ strongly dependent on the number 
of shells taken into account. 

Overall, both LDA and LDA+DMFT deliver very con¬ 
sistent results for the exchange parameters, although the 
latter are slightly smaller. Once the dynamical correla¬ 
tions are introduced, as long as the topology of the Fermi 
surface is unchanged, the main effect is to produce car¬ 
rier mass renormalization. A qualitative explanation for 
the overall decrease of the ,/ tJ ’s in LDA+DMFT can be 
found in Ref. 51, where a direct link between the total 
exchange coupling Ji and the renormalisation factor Z 
is established. However, in a multi-orbital case, differ¬ 
ent (by symmetry) orbitals are characterized by different 
renormalization factors Z and therefore there is no simple 
scaling relation between the overall exchange couplings 
extracted from LDA and LDA+DMFT. In Appendix C 
we show the computed Z-factors for each 3 d orbital cen¬ 
tered on the atom sitting either in the surface, subsurface 
or in the middle layer. One can see that in case of Fe 
slab the Z-factors for the majority and minority states 
are very different, reaching the maximal difference for 
surface electrons (being 0.53 and 0.73, respectively). A 
more detailed information about the mass enhancement 
as a measure of the strength of the many-body effects in 
all studied surfaces can be found in Appendix C. 


1. Different surface directions of bcc Fe 

In order to see whether the AFM exchange coupling 
between the two nearest neighbors at the surface of Fe 
can be observed in different surface directions, we have 
performed some additional calculations for the directions 
of (110) and (111) of bcc Fe. The obtained exchange pa¬ 
rameters between the atom at the surface and the ones 
in any layer (Jij) are shown in Fig. 6. As clear from the 
left panel of Fig. 6, the strong FM coupling happens for 
the NN atoms at the surface of (110) direction as well as 
between the second and third NNs. In contrast, such a 
strong coupling has not been observed between the near¬ 
est atoms at the surface of (111) direction. The reason 
is that in this case, the closest atoms at the surface are 
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FIG. 6. Layer-resolved exchange parameters for bcc Fe slab 
in (110) and (111) directions (top and bottom panels respec¬ 
tively). Here, we only show the results for the case when atom 
i is located at the surface (Jij) . The solid lines indicate LDA 
results while the dashed lines represent LDA+DMFT results. 
The layer numbering in the legend starts from surface denoted 
by 1, subsurface denoted by 2 and so on. 


the fourth NN of each other, too far to show significant 
coupling. As a conclusion, our calculations show that the 
AFM coupling can only be seen at the surface of (001) 
direction. 



FIG. 7. Layer-resolved projected density of states of 3 d or¬ 
bitals of Co slab for the atoms sitting at the surface, sub¬ 
surface and the innermost layer for majority and minority 
spin components in LDA (left panel) and in DMFT approach 
(right panel). 


B. Co 

Next, we have considered a slab of hep Co containing 
15 layers repeated in the (0001) direction. In Fig. 7 the 
PDOS for the innermost layer, subsurface and surface 
atoms are shown. The results obtained through both 
LDA and LDA+DMFT are reported. One can see that 
the PDOS for atoms in the innermost layer and in the 
subsurface are similar around the Fermi level, although 
their overall shapes are slightly different. For instance, 
the first peak below Ep for majority spin states has lower 
intensity and located at lower energies for atoms at the 
subsurface, in comparison with atoms in the innermost 
layer. 

Similar conclusions can be drawn for the surface. The 
majority spin PDOS for atoms at the surface shows quite 
similar features around Ep as for atoms in the innermost 
layer. This means that the strong ferromagnetic behav¬ 
ior is found both for the surface of Co like as well as the 
bulk. Nevertheless, the sharp peak below the Fermi level, 
which mostly arises from d x 2_ y 2 and d xy contributions, is 
slightly suppressed and shifted to lower energies for atoms 
at the surface. Conversely, at the surface, there is a slight 
increase in the spectral weight around the Fermi level for 
the minority spin states. This weight is mainly due to or¬ 
bitals with d yz and d xz symmetry. In addition, as we saw 
for the case of Fe, this will result in larger values for the 
local exchange field (£ s ) at the surface and consequently 
in larger exchange integrals. Finally, Fig. 7 also shows 
that the PDOS obtained from LDA and LDA+DMFT 
show similar features around the Fermi level. Hence, one 
would expect a similar trend in the asymptotic behavior 
of the exchange parameters obtained within these two 
methods. 

In Fig. 8, layer-resolved exchange parameters are dis¬ 
played for the physically most interesting layers. Results 
from both LDA and LDA+DMFT are reported. In con¬ 
trast to the Fe slab, a relatively faster decay of the ex¬ 
change parameters can be seen for the Co slab. This 
is due to that the RKKY character is less effective in 
strong magnets, e.g. as pointed out in Ref. 52. However, 
itinerant magnets, in general, are not perfect strong mag¬ 
nets, due to the hybridization between d orbitals and sp 
states. Our statement of strong ferromagnetism in this 
paper should be viewed as describing a situation when 
the majority spin states of the DOS is completely filled 
and hence pushed below the Fermi level. 

In agreement with the results for the Fe slab, we ob¬ 
tained that the interlayer exchange parameters between 
NNs are substantially smaller in the inner layers than in 
layers close to the surface (for comparison, see the green 
lines in the middle panel of Fig. 8 and the pink lines in 
the left panel). However, in contrast to the Fe surface, 
there is a strong intralayer FM coupling between the NN 
at the Co surface, both in LDA and in LDA+DMFT 
(light green lines in the right panel). An analysis of in¬ 
dividual orbital contributions to these exchange parame¬ 
ters reveals the there are strong FM contributions arising 
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FIG. 8. Layer-resolved exchange parameters (Jij) for 15-layer hep Co (0001) slab for the case when atom i is located in the 
innermost layer (left panel), in the subsurface (middle panel) and at the surface (right panel). The solid lines indicate LDA 
results while the dashed lines represent the LDA+DMFT results. The layer numbering in the legend starts from surface denoted 
by 1, subsurface denoted by 2 and so on. The innermost layer is denoted by 8. The interaction between the surface and layer 
9 and further are small and are not shown in these plots. 


from all 3 d orbitals where the d yz — d yz and d xy — d xy 
contributions are the strongest. 

Despite the strong intralayer and interlayer FM cou¬ 
pling between atoms at the surface, the associated total 
exchange interaction (Ji) is still smaller than those ob¬ 
tained for the inner layers, due to the lower coordina¬ 
tion number. As for Fe, the consequence is, qualitatively, 
in agreement with the conclusions reported by Turek et. 
al [32], but the magnitude is very dependent to the num¬ 
ber of shells included in the calculation of Jj. 

As seen in Fig. 8, both LDA and LDA+DMFT ap¬ 
proaches deliver quite similar results for the exchange 
parameters. The only difference, as we saw for the case 
of Fe, is the reduction in magnitude of the Jij ’s obtained 
within LDA+DMFT approach, while the overall behav¬ 
ior is similar and the sign is the same. 



FIG. 9. Layer-resolved projected density of states of 3 d or¬ 
bitals of Ni slab for the atoms sitting at the surface, subsurface 
and the innermost layer for majority and minority spin com¬ 
ponents in LDA (left panel) and in DMFT approach (right 
panel). 


C. Ni 

Finally, we have considered a slab of fee Ni consist¬ 
ing of 15 layers repeated in the (001) direction. In Fig. 
9 the PDOS for atoms in the innermost (middle), sub¬ 
surface and surface layers are reported, for both LDA 
and LDA+DMFT. The PDOS for an atom in the inner¬ 
most layer is similar to that of an atom at the subsur¬ 
face, specially at the vicinity of Ep. Discrepancies are 
visible at higher binding energies, where the peaks be¬ 
come narrower at the subsurface. These differences are 
small. Around the Fermi level, the PDOS of the major¬ 
ity spin for atoms at the surface is similar to that for the 
innermost layer. 

Minority spin states, instead, show some differences 
between atoms at the surface and in the innermost layer. 
These differences originate mainly from that the surface 
exhibits a larger contribution of d z 2 and d x 2_ y 2 states to 
the Ep. 

Similarly to Fe and Co, the PDOS obtained via LDA 
and LDA+DMFT are rather similar around the Fermi 
level, but possess stronger differences at higher excita¬ 
tions energies. 

Layer-resolved exchange parameters for Ni are re¬ 
ported in the top panel of Fig. 10, for both LDA and 
LDA+DMFT. Similar to the case of hep Co, a rela¬ 
tively fast decay in exchange parameters with distance 
has been observed. This is consistent with the less pro¬ 
nounced RKKY character reported for bulk strong fer- 
romagnets [52]. However, the magnitude of the coupling 
for Ni is obtained to be about three times smaller than 
those for Co. Like in the case of Fe and Co, the NN 
interlayer exchange parameters are larger for the layers 
close to the surface (for comparison see the pink lines in 
the left panel and green lines in the middle panel). As 
we have seen for the Co surface, there is a FM in-plane 
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FIG. 10. Layer-resolved exchange parameters (Jij) for 15-layer fee Ni(001) slab for the case when atom i is located in the 
innermost layer (left panel), in the subsurface (middle panel) and at the surface (right panel). The solid lines indicate LDA 
results while the dashed lines represent the LDA+DMFT results. The layer numbering in the legend starts from surface denoted 
by 1, subsurface denoted by 2 and so on. The innermost layer is denoted by 8. The interaction between the surface and layer 
9 and further are small and are not shown in these plots. 


exchange coupling between the NN at the surface of Ni 
both in LDA and LDA+DMFT approaches (green lines 
in the right panel). 

Interestingly, the sum of all exchange parameters for 
atoms at the surface (Ji) is significantly larger than the 
corresponding sum (J 2 ) at the subsurface, for both LDA 
and LDA+DMFT. This amount is 47.96 (41.87) meV in 
LDA (LDA+DMFT) for the surface versus 46.36 (41.35) 
meV for the subsurface. This might seem to be in con¬ 
trast to the fact that a lower coordination number should 
lead to a lower total exchange parameter. However, equa¬ 
tion (2) shows that the dependence on the exchange split¬ 
ting is more relevant, which explains our results for Fe, 
Co and Ni. We should also mention that here the total 
exchange parameters are evaluated inside a shell of 10 A 
radius. 


Spin and orbital moments 

As mentioned above, all calculations presented so far 
were performed in the scalar relativistic limit. It is in¬ 
teresting to analyze the influence of spin-orbit coupling 
on the magnetic properties of Fe, Co and Ni surfaces. 
Therefore, we have performed additional relativistic cal¬ 
culations including spin-orbit coupling corrections, whose 
results are reported in Table II. Our results for the in¬ 
nermost layers are in qualitative agreement with a prior 
study for bulk [28] and for the surface [17, 58, 59]. Com¬ 
paring the results of the Fe slab in Table II between LDA 
and LDA+DMFT approaches, reveals that in general 
DMFT tends to improve the results of LDA for the slabs 
both in the bulk-like and in the surface regions. This 
improvement is not only on the magnitude of spin and 
orbital moments but also in their ratio (/r \/\x s ) which for 


bcc Fe is about 0.023 (0.028) for LDA (DMFT), as can be 
deduced from Table II. Experimental value reported for 
bcc Fe bulk is about 0.037 [53]. For the surface spins, the 
enhancement of [i\/\x s is even larger (0.033 in LDA and 
0.046 in DMFT) thanks to the more pronounced orbital 
polarizations at the surface rising from more localized 
states. For hep Co and fee Ni similar conclusions can 
be drawn. DMFT enhances the value of the orbital mo¬ 
ment, both for bulk and surface atoms. For fee Ni the 
experimental values of the bulk are in good agreement 
with theory, while for Fe and Co theory underestimates 
the value of the bulk orbital moment with 15-25%. How¬ 
ever, DMFT provides a systematically better approach 
to investigate the orbital moments of these materials, at 
least judging from the bulk values. Unfortunately sur¬ 
face orbital moments are not frequently reported for these 
materials, and we list in Table II one measured value of 
fee Co (on a Cu 001 substrate) that shows enhancement 
compared to bulk values. On the other hand, spin mo¬ 
ments obtained from LDA are marginally modified by the 
dynamical correlations. 


Finally, we should mention that our results for (/Lq//z s ) 
are not quantitatively comparable to some of the recent 
experimental data based on electron magnetic circular 
dichroism (EMCD), which have reported higher values 
for this ratio (0.08+0.01 for bcc Fe in Ref. [60] and 
0.14+0.03 for hcp-Co in Ref. [61]). However, we found a 
closer agreement between theoretical results and the ex¬ 
perimental data based on X-ray magnetic circular dichro¬ 
ism (XMCD) as shown in Table [53-55]. 
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TABLE II. Layer-resolved spin (/r s ) and orbital moment (/r 0 ) for Fe, Co and Ni slabs using LDA (DMFT) approach including 
spin-orbit coupling corrections as well as the experimental values. 



calc 

M's 

Fe 

Co 

, calc 

A 4 i 

calc 

H's 

Ni 


, calc 

Ad 

calc 

H's 

, calc 

A h 

Surface 

2.92 (2.94) 

0.107 (0.122) 

1.75 (1.79) 

0.086 (0.122) 

0.75 (0.77) 

0.066 (0.080) 

Subsurface 

2.32 (2.36) 

0.057 (0.067) 

1.68 (1.73) 

0.077 (0.111) 

0.67 (0.68) 

0.055 (0.067) 

Middle 

2.20 (2.24) 

0.051 (0.063) 

1.65 (1.70) 

0.076 (0.108) 

0.63 (0.65) 

0.047 (0.057) 


calc 

H's 

calc 

Hi 

calc 

H's 

calc 

Hi 

calc 

H's 

calc 

H'l 

Bulk 

2.15 a 

0.080 a 

1.52 a 

0.140 a 

0.51 a 

0.043 a 

Bulk 

2.08 b 

0.092 b 

1.52 b 

0.147 b 

0.52 b 

0.051 b 

Bulk 

1.98 c 

0.085 c 

1.62 c 

0.154 c 

0.65 d 

0.055 d 

Bulk 



1.86 e 

0.130 e 



Bulk 



1.72 f 

0.134 f 



Surface 



1.92 f 

0.234 f 



a Ref [53]. 

b Ref [54], 

c Ref [55]. d Ref [56]. 

e Ref [57]. 

f Ref [19]. These values are for 

fee Co. 


IV. CONCLUSIONS 

In this work we have investigated the interatomic ex¬ 
change of bcc Fe, lrcp Co and fee Ni, as one comes closer 
to the surfaces from the bulk region. Our theoretical 
method is based both on the local spin-density approxi¬ 
mation, as well as dynamical mean field theory, in which 
dynamic correlations are treated explicitly. We have used 
a slab geometry in these studies, and found that for the 
central layers of the slabs, bulk like moments and ex¬ 
change parameters are found in all three studied cases. 
As one approaches the surface region from the bulk, will 
observe a general trend of enhanced spin and orbital mo¬ 
ments, both in LSDA and in DMFT. In fact, the dif¬ 
ference between results of LSDA and DMFT is rather 
minor for these systems, at least when it comes to spin 
moments and interatomic exchange parameters. For the 
orbital moments we observe somewhat larger differences 
between LSDA and DMFT results, and we find that the 
latter compare in general better to experiments, where a 
comparison is possible. 

Inspection in more detail of the interatomic exchange 
interactions reveal a general trend of enhanced values at 
the surfaces. We find that this is primarily driven by the 
increased exchange splitting of the surface states, some¬ 
thing which is caused by the reduced coordination num¬ 
ber of surface atoms. Hence, the experimental observa¬ 
tion of lower ordering temperatures of surfaces, which is 
a general phenomenon, is not caused by a reduction of 
the interatomic exchange interactions of the surfaces. In 
contrast the surface exchange interactions are enhanced. 
However, when coupled to an effective spin-Hamiltonian, 
of Heisenberg type or similar, the reduced coordination of 
surface atoms reduced the local Weiss field of the surface 
atoms, which makes them more susceptible to thermal 
fluctuations. 

Finally we have analyzed symmetry resolved aspects 
of nearest neighbor interactions of surface atoms of bcc 


Fe, and found that some of these interactions are ferro¬ 
magnetic whereas some are antiferromagnetic, and that 
summed over all symmetry components, the nearest 
neighbor exchange interaction of surface atoms is anti¬ 
ferromagnetic. The magnetic order of the Fe surface is 
nevertheless ferromagnetic, due to strong ferromagnetic 
coupling to subsurface atoms. We argue however that 
the antiferromagnetic surface interactions of bcc Fe are 
inherent, and should be an avenue to tune complex mag¬ 
netic structures of mono-atomic overlayers of Fe on bcc 
substrates. 
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Appendix A: The effect of U on J ij 

At the surface, because of the less effective screen¬ 
ing, the value of the Hubbard U is expected to increase. 
In this regard, we performed series of calculations using 
higher values for the surface and subsurface atoms to see 
the impact of U on the exchange parameters. For Fe, we 
have used U = 3 eV for atoms at the surface and U = 2.8 
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eV for those in the subsurface, while value for the inner 
layers is kept fixed to 2.3 eV. The outcome of these calcu¬ 
lations, together with the results obtained with a uniform 
U value (2.3 eV) for all atoms, are shown in Fig. A.l. It 
is evident that larger U values for surface atoms result 
in a small uniform reduction of the parameters Jij , but 
the trends are unchanged (comparison of the dashed lines 
with the solid lines of Fig. A.l). From the extent of these 
changes, we conclude that the overall behavior and more 
in particular the sign of the coupling do not change if U is 
varied within a reasonable range. Finally, similar conclu¬ 
sions can be obtained from analogous calculations for Co 
and Ni, which confirms that the values of the exchange 
parameters are rather robust with respect to the choice 
of the U value. 




FIG. B.l. Layer-resolved exchange parameters (Jij) for the 
surface atom of 7-layer Fe slab (atom i is located at layer 1). 
The solid (dashed) lines indicate the results without (with) 
updating the electron density. The layer numbering starts 
from surface denoted by 1, subsurface denoted by 2 and so 
on. 

the transition metal surfaces. This may be important 
for future investigations of thin films deposited on a sub¬ 
strate, where computational efficiency is going to be of 
primary importance. 

Appendix C 


FIG. A.l. Layer-resolved exchange parameters (Jij) for the 
surface atom of bcc Fe (atom i is located at layer 1). The solid 
lines indicate the results obtained for U = 2.3 eV for all atoms, 
the dashed lines represent the results for layer-dependent U 
values (see the text). The layer numbering starts from surface 
denoted by 1, subsurface is denoted by 2 and so on. 


Appendix B: Charge self-consistency 

The results presented in the main text refer to cal¬ 
culations where the electron density is kept fixed to its 
LDA value, and the local correlation effects affect the 
results only through the self-energy function. We how¬ 
ever performed several calculations to analyze the role of 
complete self-consistency over the electron density. We 
illustrate these results for a 7-layer slab of Fe(001). The 
obtained exchange parameters for one atom at the surface 
are reported in Fig. B.l, with and without updating the 
electron density. One can observe only small variations 
in the absolute magnitude of the parameters . These 
differences amount to only a few percent. Similar calcula¬ 
tions have been performed for 7-layer slabs of Co and Ni, 
and lead to similar results. Thus, we conclude that the 
effects of charge self-consistence within the LDA+DMFT 
scheme are negligible for treating magnetic properties of 


In this appendix we show the results for the calculated 
orbital-resolved renormalisation factors Z in slabs of Fe, 
Co and Ni. Z ^ denotes the inverse of an effective mass 
enhancement for the correlated orbital m and is defined 
as: 



du> 



(Cl) 


where S^ m (w) is the self-energy projected on the orbital 
m with the spin projection a, u> = 0 corresponds to Ep. 
In general, Z-factors are good measures of the strength 
of the correlation effects. To have a compact description 
of the latter, we have also calculated the average Z% 
per spin channel. Z° vg was computed using the following 
expression: 


avg 


£m Z mN m (Ep) 

E m N m (E F ) ’ 


nii — d z 2 ,..., d X y (C2) 


where N^(Ep) denotes the partial density of state at 
the Fermi level of a particular state. Thus, the average 
Z-factor is a weighted sum of orbital-resolved renormali¬ 
sation factors. The weight of each orbital is defined by its 
relative contribution to the spectral weight at the Fermi 
level. 

Calculated values of orbital-resolved Z-factors as well 
as their average values are shown in Table C.l. An in- 
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TABLE C.l. Orbital-resolved and average (avg.) renormalization factors Z for Fe, Co and Ni slabs. 


i Fe | 


Majority Spin 

Minority Spin 


d z 2 

d x 2_ y 2 

dy Z 

dxz 

dxy 

avr 

d z 2 

d x 2_ y 2 

dy Z 

d xz 

dxy 

avr 

Surface 

0.57 

0.54 

0.51 

0.51 

0.54 

0.53 

0.76 

0.72 

0.74 

0.74 

0.72 

0.73 

Subsurface 

0.73 

0.72 

0.72 

0.72 

0.72 

0.72 

0.78 

0.77 

0.81 

0.81 

0.81 

0.80 

Middle 

0.72 

0.72 

0.71 

0.71 

0.71 

0.71 

0.76 

0.76 

0.76 

0.76 

0.76 

0.76 


1 Co ! 


Majority Spin 

Minority Spin 


d z 2 

d x 2_ y 2 

dy Z 

dxz 

dxy 

avr 

d z 2 

d x 2_ y 2 

dy Z 

dxz 

dxy 

avr 

Surface 

0.73 

0.75 

0.74 

0.74 

0.75 

0.74 

0.73 

0.78 

0.75 

0.75 

0.78 

0.76 

Subsurface 

0.78 

0.77 

0.77 

0.77 

0.77 

0.77 

0.82 

0.81 

0.82 

0.82 

0.81 

0.81 

Middle 

0.78 

0.77 

0.76 

0.76 

0.77 

0.77 

0.82 

0.81 

0.81 

0.81 

0.81 

0.81 


1 Ni i 


Majority Spin 

Minority Spin 


d z 2 

d x 2_y2 

dy Z 

dxz 

dxy 

avr 

d z 2 

d x 2_ y 2 

dy Z 

dxz 

dxy 

avr 

Surface 

0.78 

0.77 

0.79 

0.79 

0.78 

0.78 

0.69 

0.69 

0.68 

0.68 

0.71 

0.69 

Subsurface 

0.82 

0.82 

0.80 

0.80 

0.81 

0.81 

0.77 

0.75 

0.75 

0.75 

0.75 

0.75 

Middle 

0.81 

0.81 

0.80 

0.80 

0.80 

0.81 

0.76 

0.76 

0.74 

0.74 

0.74 

0.75 


spection of the results suggests that overall the correla¬ 
tion effects in all studied systems are not very strong (as 
a limiting case, Z= 1 corresponds to LDA). In Co and 
Ni slabs the renormalisation factors for all 3 d orbitals 
are similar and their values lie in the range between 0.7 
and 0.8. It is also seen that the many-body effects are 
the most pronounced for surface electrons. The largest 
renormalisation effects are found for Fe slab. In partic¬ 
ular, majority-spin electrons of the surface atoms expe¬ 


rience an almost twice mass enhancement. Their overall 
spectral weight at the Ep is relatively small (see. Fig. 
3) and thus these quasiparticles are more sensitive to the 
addition of the self-energy. 

Here we emphasize again the fact the differences in the 
^-factors for majority and minority electrons are quite 
substantial. However, within a particular spin channel 
the orbital-resolved Z-factors show relatively smaller de¬ 
viations from the average value. 
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